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Abstract 

We have evaluated the nonmesonic weak decay spectra within the independent-particle shell- 
model, and compared them with the recent measurements of: i) the single and double coincidence 
nucleon spectra in ^^C performed at KEK, and ii) proton kinetic energy spectra in AHe, ^Li, 
ABe, A B, A^C, "C, a'^N and i^O done by FINUDA. Based on this comparison we argue that 
the extraction from the data of the three-body ANN — i- nNN induced decay rate, as done in 
these works, could be questionable. 
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1. Introduction 

The weak decay rate of a A hypernuclcus can be cast as [1] 

Tw +^NM, (1) 

where Tm is decay rate for the mcsonic (M) decay A ttN, and Tnm is the rate for 
the nonmesonic (NM) decay, which can be induced either by one bound nucleon (lA'^), 
ri(A7V nN), or by two bound nucleons (27V), r2(AA^A^ nNN), where N = p,n 
i.e., 

Ttva/ = Ti + Ti = Tp + r„, T2 = r„„ -(- Tnp + Tpp. (2) 
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While the M and liV-NM decays have been seen experimentally in the pioneering 
measurement performed more than 50 years ago by Schneps, Fry, and Swami [2], the 
experimental observation of the 2 A^-NM decay, which was predicted by Alberico et al. [3] 
in 1991 (sec also Ref. [4]), has been reported only quite recently at KEK [5], and at 
FINUDA [6]. Both groups announced a branching ratio T2/rNM '--^ 25 — 30 %. The first 
group obtained this result from the single and double coincidence nucleon spectra in ]^C, 
and the second one from proton kinetic energy spectra in ^He, ^Li, ^^e, ]^C, ]^C, 
j^N and )^0. A branching ratio for the 2A^-NM decaying channel of such a magnitude 
is consistent with the prediction made by Bauer and Garbarino [7], but it is very large 
in comparison with the upper limit r2/rvK < 0.097 (95% CL), established previously 
from the single-particle proton and neutron kinetic energy spectra in \iie at BNL [8]. In 
this experiment it was also found that the effect of the final state interactions (FSIs) is 
relatively small, its upper limit being Tf^jj/Tw < 0.11 (95% CL). ^ 

Quite recently, Bauer and Garbarino [10] have obtained good agreement with KEK 
data [5] , considering both the one- and the two- nucleon induced decays in the framework 
of the Fermi Gas Model (FGM). These authors have also analyzed the proton spec- 
trum measured at FINUDA [6], but no theoretical study of the remaining spectra has 
been done so far. The aim of the present work is twofold. First, we confront the results 
of the independent particle shell model (IPSM) for the lA^-NM decay with the above 
mentioned experiments. Second, based on this comparison we expose our doubts about 
the reliability of the procedures followed in these works to extract the ratio r2/rArA/ 
from data. 
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Fig. 1. Schematic representation of a) one- nucleon, and b) two- nucleon induced decays in A-hypernuclei 
when described by the interplay of weak {W) and strong (S) interactions through the exchanges of 
nonstrange-mesons 7r,p, cj, and rj, and strange-mesons K, and K*. 

The schematic representation of the two channels, when the decay dynamics is ac- 
counted for by the one meson-exchange (OME), is shown in Fig. 1. This is the most fre- 
quently used model for handling the NM-decay, and includes the exchanges of nonstrange- 
mesons 7r,p, w, and 77, and strange-mesons K, and K*. It is based on the original idea 
of Yukawa that the NN interaction at long distance is due to the one-pion-exchange 
(OPE), the dominant role being played by the exchange of pion and kaon mesons. 



^ In fact, the ^Hc spectra arc accounted for reasonably well theoretically by considering the lA^-NM 
decay mode only [9], 
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The OPE potential was verified quantitatively by the Nijmegen partial wave analysis of 
the NN scattering in the elastic region [11]. i.e., at distances larger than the minimal de 
Broglie wavelength 1/ y/rriT^yi ^ 0.5 fm corresponding to the pion production threshold. 
The verification of other meson exchanges is less straightforward, and the uncertainties 
in the baryon-baryon-meson (BBM) coupling constants could be quite sizeable since they 
are not constrained by experiments. To derive them in the strong sector (S vertices in 
Fig. 1) the SU{3)f (flavor) symmetry is utilized. In the weak sector (W vertex in Fig. 
1) the BBM parity- violating couplings are obtained from the SU{6)w (weak) symmetry, 
while the parity conserving ones are derived from a pole model with only baryon pole 
resonances [12]. 

The uncertainties regarding the BBM coupling constants are amplified still further 
by the short-range correlations (SRCs) between the emitted nucleons nN, and nNN. 
Parreiio, and Ramos have shown that they can diminish the value of Fi by more than 
a factor of two [13]. Nothing has been said so far regarding the effect of the SRCs 
on the 2A^-NM decay. The theoretical scene becomes still more complex when effects of 
quark degrees of freedom [14,15], the 27r-exchanges [16-18], and the axial- vector ai-meson 
exchange [18] are considered. 

To derive the ratio T2/TNA1 from data, both KEK [5] and FINUDA [6] teams have 
included the FSIs in their analysis. The first group has used them to gauge their exper- 
imental data on single and double coincidence nucleon spectra in ]^C. For this purpose 
they have used the Intranuclear Cascade (INC) code, which was developed by Ramos 
et al. [19] to follow the fate of nucleons produced by the primary lA^-NM and 2A^-NM 
decays. On the other hand, the FINUDA's group [6], in their analysis of proton kinetic 
energy spectra in several hypernuclei, have assumed that the basic role of the FSIs is 
to shift the transition strength from high energy to low energy, without substantialy 
modifying the total number of protons. 

The above mentioned experiments, together with several others performed during the 
last decade [20-28], represent very important advances in our knowledge about the NM 
decay. Explicitly, they are: 1) new high quality measurements of single-nucleon spectra 
Sn{En), as a function of one-nucleon energy E^, and 2) first measurements of the two- 
particle-coincidence spectra, as a function of: i) the sum of kinetic energies En + En = 
EnN, SNiEnN), h) the opening angle 0,ijv, SNicosOnN), and iii) the center of mass (cm.) 
momentum P„n = \PnN\-, SjsfiPnN) with FnN = Pn + Pat- On the theoretical side this 
implies a new challenge for nuclear models which now have to explain, not only the lA'^- 
and 2Af-NM decay rates, but also the shapes and magnitudes of all these spectra, testing 
in this way both the kinematics and the dynamics. 

The measured spectra are obtained by counting the numbers of emitted nucleons ANjv 
within the energy bin AE = 10 MeV, or the angular bin Acos0 = 0.05 always corrected 
by the detection efficiency. Here we take advantage of the fact that in the KEK exper- 
iments [5,25,29] ANjv are normalized to the number of NM processes ^nm, while for 
the FINUDA proton spectra [6,30] we have at our disposal also the number of produced 
hypernuclei N^. Therefore, we can exploit the following relationships 

AN^ ^ AT^ AN^ ^ AT^ 
NjvA/ Fata/ Nvk Tw 

where AF^^ is the emission rate of protons within the experimentally fixed bin. Thus, 
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where N^r = X^I^i ^N^r is the total number of NM events decaying to the mode N, and 
= J2^i ^r^^, is the corresponding decay rate; m is the number of bins. We note 
that the last relation in (4) agrees with Eq. (13) in Ref. [8], and that accounts for 
both the lA^-NM and 2A^-NM decays as well as for the FSIs. 
The correspondence between theory and data is 
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For is Tnm = 0.95 ± 0.04 [5], and as it was pointed out in Ref. [28] 



Tw{A) = (0.990 ± 0.094) + (0.018 ± 0.010) A, 



(6) 



for all measured hypernuclei in the mass range A = 4 — 12. The theoretical decay widths 
are AF^ = Sn{En)AE, etc. Below we briefly sketch the corresponding expressions for 
spectral densities Sn within the IPSM. 



2. IPSM for the lAT-NM decay 

Within the IPSM [9,31-40]: i) the initial hypernuclear state is taken as a hyperon A in 
single-particle state j'a = lsi/2 weakly coupled to an {A — 1) nuclear core of spin Jc, i.e., 
I J/) = \{JcjA)Ji), ii) the nucleon {N = p, n) inducing the decay is in the single-particle 
state Jn (j = nlj), iii) the final residual nucleus states are: \ Jf) = \{JcjN^)JF), iv) the 
liberated energy is 

A^^ = A + £A + e]^, (7) 

where e's are single-particle energies, and A ~ Ma — M, and v) the cm. momenta, and 
relative momenta of the emitted particles are: 



A 



PnN ^^{A- 2)(2MA^^ ~pI~ pI), p^n - JMA^^ - ^^^l^f nN- (8) 



^(A-2){2MA],~pI~p%), PnN^sj- 
It follows that the lA^-NM decay rate reads 

FAr=^F]^; F]^ J lj^{pnN,PnN)d^lnN, (9) 
3 

where (if2„jv is the phase space factor, which depends on the spectra that one is interested 
in, and 

J=j + l/2 

lUpnN,PnN)= ^ F^,, ^ T^^^ (p„Jv)Oi (P„iv). (10) 

J=\]-l/2\ L 

The information on nuclear structure is contained in the spectroscopic factor 
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^^j = ^/"'Ek^^ii(«L«]a), 11^^)1'- (11) 



where J = -\/2 J + 1. The values for J/, and Jc are taken from experimental data and 
for hypernuclei of interest here are listed in Tables I of Ref. [37]. The resulting factors 
are presented in Tables II of the same paper. 
The kinematics arc enclosed in: i) the phase space factor dQnN, and ii) the overlap 

OLiPnN)^ J R^dRjLiPnNR)R0Lib/V2,R), (12) 

between the cm. radial wave functions Rql of the bound particle and of the outgoing 
particle, where b is the harmonic oscillator size parameter. 
The decay dynamics is contained in 

UjLiPnN) - J2 \^(PnN] ILXS JT ] J Aj N JtAN)\^ , (13) 
SIXT 

where X = I + L, with I and L being, respectively, relative and cm. angular momenta, 
T = {TMt,Mt = mt^+mtr^}, and tAAr = {^a = 1/2, mt^ = -1/2, = 1/2, mt„}, with 
mt = 1/2, and rrit^ = —1/2, where we have assumed that the AN — >■ nN interaction 
occurs with the isospin change AT =1/2. Moreover 

M{pnN;lLXSJT;jAjNJtAN) = ^ [1 - (-)'+''+^] 

X {lLXSJT\V{pnN)\jAjNJtAN), (14) 

where (and henceforth) the ket |), unlike |), indicates that the state is not antisym- 
metrized, and V{pnN) is the transition potential. 

To find the kinetic energy and opening angle d„N spectra it is convenient to express 
the single-particle transition rate as 



■ N 



(A- 2)— — ^ / dcOs6'„Ar / dENd ^ EnIjj^{pnN,PnN), (15) 



with 

E'j^ = iA^ 2) (A - 1)A,„ - En[{A - 1)2 - cos^ 0„^], (16) 



and 
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(18) 



Throughout the integration one has to enforce the condition E'j^ > En cos^ ■ It 
might be worth noticing that, while E'j^ does not have a direct physical nicaning, En is 
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the energy of the neutron that is the decay-partner of the nucleon N with energy E^. 
The corresponding spectra read 



+1 , 

dV SM"^ f IE 

Sn{En)^-T^^ {A~2) V / dcOSenNj EnI,^{pnN,PnN): (19) 



In _i 

and 



SNicOSe^N)^:: ^ = (A-2) V / dENJ-^E„Ij,{p„N,PnN). (20) 

As one-proton (one-neutron) induced decay prompts the emission of an np (nn) pair the 
total neutron kinetic energy spectrum is: 

S,itiEn) = , "* = Sp{En) + 2Sn{En), (21) 

where T„t — Tp + 2r„ is the total neutron transition rate. 
Similarly, we get 

SnIEun) = v/^(A-2)3 V J (A,„ - EnN){EnN ~ Aj„)I^(p„Ar, P„jv), (22) 

TT — ' ' 

Jn 

for the kinetic energy sum EnN = En + E^, and 



5iV (PnTV ) = ^ y ^ J2 PnN ^ Pf. - Kn'^,. (PnN . Pun) . (23) 

jN 

for the cm. momentum PnN spectra, with 



A~2 , „ „ A-2 



Aj„=A,„^-, and P,„ =2y^-MA,„ (24) 

being, respectively, the maximum values of EnN^ and PnN for each jn- 

The outline of the numerical calculation is the following: a) The transition potential 
for the emission of the nN pair, V{pnN) in Eq. (14), is described by two OME models 
with the weak coupling constants from Ref. [12,13], namely: Ml) only the tt + A' exchange 
is considered, and M2) the full 'K-\-K + r] + p-\-uj + K* is taken into account, b) The 
initial and final SRCs, as well as the finite nucleon size effects are included in the same 
way as in our previous works [9,31-38], and c) The parameter b in Eq. (12) is evaluated 
as in Ref. [32], i.e., b = l/^/h^JM^, with hu ^ A5A-^/^ - 25^-2/3 ^j^y. 



3. KEK Experiment 



In Figs. 2, 3, and 4 we compare the KEK [5,29] experimental data with the IPSM 
calculations. As seen from Fig. 2 the calculated proton and neutron kinetic energy spectra 
agree fairly well with data for En > 40 MeV, and specially when the it + K potential is 
used, in spite of not considering the FSIs. This is not the case, however, for the angular 
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Fig. 2. Comparison between tlie experimental and theoretical kinetic energy spectra for neutrons (lower 
panel), protons (middle panel), and the sum of both (upper panel). The KEK experimental data are 
from [5,29], and the relation (5) has been used. The theoretical results have been evaluated from the 
relations Arp(E) = Sp{E)AE, and Ar„t{E) = [Sp{E) + 2S„(E)]AE, where AE = 10 MeV is the 
experimental energy bin, and Sn{En) are given by (19). Two OME potentials have been used: Ml) only 
IT + K exchange is considered, and M2) full n + K + ri + p + Lj + K* exchange is taken into account. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 

COS 

Fig. 3. Comparison between experimental and calculated opening angle correlations for neutron-neutron 
pairs (lower panel), proton-neutron pairs (middle panel), and the sum of both (upper panel). The KEK 
experimental data are from [5,29], and the relations (5) have been used. The theoretical results have 
been evaluated from the relations Arjv(cos 6„jv) = Spf{cos9npf)Acos9, where Acos(9 = 0.05 is the 
experimental opening angle bin, and Sm {cos 6„n) are given by (20). The same OME potentials have 
been used as in Fig. 2. 
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distributions, and the cm. momenta of the nN pairs shown in Figs. 3, and 4, where the 
theory overestimates the data quite signifieantly, and particularly for the np pairs. It 
could be worth noticing that both potentials yield very similar results for the nn pairs. 




100 200 300 400 500 600 
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Fig. 4. Comparison between experimental and calculated cm. momentum for neutron-neutron pairs 
(lower panel), proton- neutron pairs (middle panel), and the sum of both (upper panel). The KEK 
experimental data are from [5,29], and relations (5) have been used. The theoretical results have been 
evaluated from the relations Arjvl^ViJv) = Spf{Pn]sf)^P, where AP = 50 McV/c is the experimental 
center of mass momentum bin, and SAr(PnJv) a^re given by (23). The same OME potentials have been 
used as in Fig. 2. 

The experimental decay rates Tn are obtained by summing over all AFjv in Figs. 2, 
3, and 4. The results are shown in the first three rows of Table 1. In the following two 
rows the results of the IPSM calculation are exhibited. No energy or angular cutoffs are 
done on the calculated results, since their effects are insignificant. It is self evident that 
the theory overestimates the data and in particular so in the case Ml. We are aware that 
a more elaborate description of the SRCs, as performed in Ref. [13], could improve the 
agreement between data and theory for transition rates which are shown in Table 1. In 
this regard the FSIs could also help, which in addition are able to modify the spectra by 
redistributing the calculated transition strengths. On the other hand the new 2iV-NM 
channel can only increase the theoretical results but never decrease them, and thus the 
discrepancy with experiment would be still more pronounced. 

The measurement of F2 in Ref. [5] is based on INC-liV calculations that overestimate 
the data above thresholds En > 30 MeV for both single and coincidence spectra. The 
authors denominate this difference as "quenching" of data in relation to theory [5,26]. 
Moreover, they attribute it to the lack of the 2A^-NM channel, which when incorporated, 
in proportion of T2/T]yM = 0.29 ± 0.13, decreases the calculated decay rate TjyM and 
yields agreement with data. Several observations are pertinent here: 

1) Usually the word "quenching" has a different meaning in nuclear physics. The 
emblematic example is that of Gamow- Teller (GT) strength, the experimental value 
of which, either from nuclear /3-decay or from charge-exchange reactions, is in general 
overestimated by calculations. This discrepancy is often resolved by including additional 
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Tabic 1 

Comparison between the experimental transition rates derived from Figs. 2, 3, and 4, and the IPSM 
calculation, making use of two OME potentials: Ml) exchange of tt + K only, and M2) full n + K + r} + 
p + w + K* . The results for Fig. 2 correspond to the proton and neutron thresholds Ejv > 35 MeV. 



Source 


r^i, Tp -|- r^^ 




Fig. 2 


0.419 ±0.013 


0.823 ± 0.016 1.242 ± 0.021 


Fig. 3 


0.197 ± 0.024 0.147 ± 0.022 0.345 ± 0.033 




Fig. 4 


0.197 ± 0.023 0.136 ± 0.022 0.338 ± 0.033 




Ml 


0.627 0.201 0.828 


1.455 1.656 


M2 


0.792 0.205 0.997 


1.789 1.994 



degrees of freedom; neither the total theoretical value of the GT operator ctt is modified 
nor opening of new decay or reaction channels is required. 

2) The present calculation also overestimates the data. Nevertheless, we do not dare to 
state that the NM transition operator V{pnN) is quenched, since as pointed out previously 
its evaluation involves many uncertainties, which is not the case with the GT operator. 

3) The INC sequence in the calculations of the FSIs performed so far [39,40,44,45] is 
always triggered by the primary nN or nNN nucleons. These nucleons are produced by 
the same OME dynamics as we have used here. In the KEK calculation the primary NM 
weak dynamics is totally ignored [5]. 

4) The Monte Carlo study of \^C performed in Refs. [39,40], where the production 
of primary nucleons is described according to IPSM, does not show any evidence of 
the "quenching" reported in Ref. [5]. This fact suggests that the uncertainties in the 
FSIs calculation must be carefully considered before extracting definite conclusions from 
comparison of the INC results with the experimental data. 

5) The above mentioned uncertainties embedded in the KEK description of the primary 
NM decay dynamics cannot but be augmented when the FSIs are considered. Therefore, 
the ambiguities involved in INC-liV calculation are probably as large or still larger than 
the observed " quenching" . 

Briefly, in our opinion the reckoning followed by the KEK group is not sufficiently 
robust to support the statement that the NM niViV-decay mode has been unequivocally 
observed experimentally. 

4. FINUDA Experiment 

The comparison between FINUDA [6,30] experimental data and the IPSM calculation 
for the proton kinetic energy spectra is shown in Fig. 5. One sees that the theory over- 
estimates the data significantly, particularly in the case of ^Li. Similarly as in the case 
of the KEK data, discrepancies between data and theory, are manifested also in the case 
of FINUDA transition rates Tp, as can be seen from Table 2. 

The determination of the branching ratio T2/TNM in Ref. [6] is based on: i) separating 
the total number of detected protons Np into low, and high energy regions populated by 
Np , and N^ protons, respectively, and ii) assuming that the lA^-NM decay rates of the 
primary protons N^ = N(Ap np) arc equal in these two regions, i.e., P^'^ = T^'^ — 
Pp/2, where the superscript indicates that this decay rate, unlike those in Tables 1 
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Fig. 5. Experimental data [6,30] for proton kinetic energy spectra arc compared with the IPSM results. 
Theoretical results have been evaluated from the relation ATp{E) = Sp{E)AE, where AE = 10 MeV is 
the experimental energy bin. The same OME potentials have been used as in Fig. 2. 

Table 2 

In the first three columns the comparison is done between the experimental transition rates Fp derived 
from Figs. 5 and the IPSM calculation, making use of two OMEP: Ml) full ir + K + ri + p + Ld + K*, 
and M2) exchange of tt + A" only. In the last two columns are listed the mean values of the Gaussian fit 
as obtained in Ref. [6], and calculated maxima within the IPSM, in units of MeV. 



Hypernucleus 


pFINUDA pMl 


pM2 
P 


gFINUDA gIPSM 




0.284 


0.466 


0.360 


68.5 ± 


4.1 


76.7 


iLi 


0.123 


0.531 


0.415 


76.7 ± 


5.2 


72.3 


iBc 


0.269 


0.627 


0.494 


78.2 ± 


6.2 


66.2 




0.380 


0.667 


0.527 


75.1 ± 


5.0 


68.9 


A ^ 


0.358 


0.792 


0.627 


80.2 ± 


2.1 


67.5 


fc 


0.437 


0.776 


0.614 


83.9 ± 


12.8 


67.6 




0.402 


0.821 


0.651 


88.1 ± 


6.2 


61.6 




0.320 


0.906 


0.718 


93.1 ± 


6.2 


62.0 



and 2, doesn't contain contributions coming from the FSIs and 2N-'NM decay. Moreover, 
they write 

N< _ o.5 + r2/r° + jv|-^VN° 

N 1 + r2/ro + TV^^^VNO ' ^ ' 

where N^'^-^ = N^'^^ + N^^' . The partition energies E are fixed as being the mean 
values of Gaussian-function fits of each proton spectrum from 80 MeV onwards, and the 
experimental values of the ratios R = /Np found in this way are approximated by a 
linear function of the mass number A, i.e., 
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Table 3 

Results for the parameters a and b, and the corresponding ratios V2/T'^, and r2/rNiv[; A) Ref. [6], 
B) data from Ref. [6], and the peak energies £;IPSM ffoj^ i^st column in Table 2 instead of £;FINUDA^ 
and C) same as in B but performing the fit with Eq. (30) instead of Eq. (26). 



case 


a 


h 




r2/rNM 


A 


0.654 ±0.138 


0.009 ±0.013 


0.43 ±0.25 


0.24 ±0.10 


B 


0.674 ±0.108 


-0.008 ± 0.009 


0.53 ±0.51 


0.26 ±0.10 


C 


0.135 ±0.408 


0.207 ±0.202 


-0.42 ±0.14 


-0.40 ±0.10 



R{A) = a + hA, where a = 

does not depend on A, and yields 
r2 a -0.5 V2 



o.5 + r2/r° 

l + Ta/ro ' 



0.5 



p 



1 



r 



NM 



(i-a)ro/r" + o.5- 



(26) 



(27) 



Finally, a fit for R{A) was done obtaining the values of a, and 6 shown in the row 
A of Table 3, together with the resulting T2/T^, and r2/rjvA/ derived from (27) with 
r^/Fp = 0.48 ± 0.08, which is the mean value of recent experimental results [26]. 

The manner of deducing the 2A^-NM decay rate as performed in Ref. [6] can be 
challenged from several aspects which we point out below. 



=1/2 O 
P3/2O 

1p„,0 
totafO 
total He 




20 40 60 80 too 120 140 160 

E(MeV) 

Fig. 6. Calculated ^He and ^^O spectra within the IPSM employing the t: -\- K transition potential. 
Contributions of different partial waves Si/2) P3/21 s-nd P1/2 to the total spectra are also displayed. 

1) How to arrive from (25) to (26) is not a trivial issue. One possibility is to neglect the 
last term in the denominator of (25), arguing, as was done in Ref. [6], that the FSIs tend 
to remove protons from the high energy part of the spectrum (N^'^-'^ < 0) while filling 
the low energy region (N^-^^ > 0), with the net result that 'N^^^ = N^-^^ + N^-^^ = 0. 
Therefore 

^FSI ^FSI 

R^ci+ gga+^, (28) 



(1 + F2/F0)N0 
which, when compared with (26), yields 



FSI 



bAKL 



(29) 
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Kinetic Energy (MeV) Kinetic Energy (MeV) 

Fig. 7. Experimental data and Gaussian fits from Rcf [6] for proton kinetic energy spectra, arc compared 
with the IPSM results normalized for proton energy Ep > 40 MeV. 

At first glance this sounds reasonable because the effect of FSIs should increase with A. 
But, since b — 0.009 ±0.013 one finds that the number of particles N;^"^^ per nucleon that 
are dislocated from high to low energy is less than 1% of the total number of primary 
protons produced by the 17V-NM decay, which is unrealistic. 

2) The assumption r^< = r°> = T°/2 is badly supported by the IPSM, where T°> 
becomes progressively larger than F^^ as the mass number is increased. Why it is so 
can be understood from inspection of Fig. 6, where the spectra of ^He and are 
displayed. In the case of ^^O the partial wave contributions are also presented. One sees 
that, while in ^He only the orbital Si/2 contributes, being Fp< appreciably larger than 
Fp> (Fp< = 0.55F°), for other hypernuclei also the orbitals P3/2, and pi/2 contribute. 
Their spectra are wider and localized at higher energies than the corresponding Si/2 
strength. This causes the total transition strength to be gradually shifted towards higher 
energies as A increases, becoming F^^ = 0.43Fp for ]^0. 




1.9 2 21 2,2 2.3 2.4 2.5 2.6 2 4 6 8 10 12 14 16 18 



A/Z Mass Number A 

Fig. 8. Ratio /? as a function of A (left panel), and of A/Z (right panel). In both cases the energies 
^IPSM^ and the corresponding R = N</N are used. 

3) Except for ^He, the calculated proton spectra deviate significantly from the Gaus- 
sian fits adopted in Ref. [6] for the liV-NM decay. Although bell shaped they differ both 
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in energies E and in FWHM widths. This can be observed from Table 2, and Fig. 7, where 
the IPSM results are normalized to the experimental data for proton energies Ep > 40 
MeV. Both OMEP models yield the same normalized spectra, since the spectrum shapes 
depend only on kinematics. 

4) Large differences between £;finuda^ g^^^^ ^ipsm^ shown in Table 2, give rise to large 
differences between values of ratio R in Ref. [6, Fig. 2], and our values exhibited in Fig. 
8. On the left and right panels of this figure are shown the linear fits of R given by 
(26), and by 

R{A,Z)^a + bA/Z, (30) 

respectively. The later is also a possible parametrization for (25), as suggested by expres- 
sion (28) since should be proportional to the number of protons Z [38]. The resulting 
values of a, b, r2/rp, and r2/rNM are listed listed in rows B and C of Table 3. One sees 
that a, and therefore r2/rp, and r2/rNM are very similar in cases A and B. However, 
the value of parameter b is negative in case B, which makes its physical interpretation 
still more embarrassing than in case A, since from (29) it follows that the number N;^'^^ 
becomes negative, and therefore the FSIs now move particles from low to high energy. 
Finally, in case C the ratios r2/r°, and r2/r NM turn out to be large and negative, which 
obviously doesn't make sense. 

We conclude therefore, that similarly to the KEK experiment [5], the argumentation 
followed by the FINUDA group [6] can yield several quite different results for the ratio 

r2/rjvj\/- 

5. Summary and Conclusions 

We have discussed the new KEK [5] and the FINUDA [6] experiments on the non- 
mesonic weak decay of hypernuclei in the framework of the IPSM, with the dynamics 
described by the OPE mechanism, and using two different models: Ml) with tt+K mesons 
only, and M2) that comprises the exchange of the complete pseudoscalar and vector me- 
son octets {■K^'q.K, p,uj,K*). We prefer the first one since the coupling constants NNtt, 
and A^Att are well established, and the only parameters that are still uncertain in this case 
arc the coupling strengths NNK . and NAK. Of course, except for reasons of simplicity, 
there is no justification for not including the remaining mesons. 

We found that the theory overestimates the data to a great extent. There could be 
several reasons for this, such as: i) the absence of the FSIs in the employed model, ii) 
the uncertainties inherent in the OME model, iii) the approximation used for the SRCs, 
etc. 

We argue that at present there exist many uncertainties, on both sides experimental 
and theoretical, which prevent us to draw definite conclusions regarding the 2A^-NM 
decay rate r2. In view of this situation we are developing the IPSM formalism for this 
decay channel [46] that will be confronted with the FGM, which is the only approach 
used so far [9] . Once this is accomplished it will also enable us to study the effects of FSIs 
on the threc-nucleon emission. New developments in this field of experimentation will be 
very welcome. Needless to say that triple (p, n, n), and {p,p, n) coincidence measurements 
will be extremely useful for direct measurement of r2. 
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After completing the present work we have learnt that a new derivation of the ra- 
tio r2/rNM has been done at FINUDA. based on the analysis of the [tt^ ,p,n) triple 
coincidence events [47] . The result is similar to the previous one [6] . 
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